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^^ , Abstract. We show that the asymptotic iteration method converges and yields 



accurate energies for a perturbed Coulomb model. We also discuss alternative 



C I perturbation approaches to that model. 
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1. Introduction 

The asymptotic iteration method (AIM) is an iterative algorithm for the solution of 
Sturm-Liouville equations [1,2]. Although this method does not seem to be better than 
other existing approaches, it has been applied to quantum-mechanical [3-5] as well as 
mathematical problems [6]. For example, the AIM has proved suitable for obtaining 
both accurate approximate and exact eigenvalues [1-6] and it has also been applied to 
the calculation of Rayleigh-Schrodinger perturbation coefficients [5,7]. 

Recently, Barakat applied the AIM to a Coulomb potential with a radial polynomial 
perturbation [5] . By means of a well-known transformation he converted the perturbed 
Coulomb problem into an anharmonic oscillator. Since straightforward application of 
the AIM exhibited considerable oscillations and did not appear to converge Barakat 
resorted to perturbation theory in order to obtain acceptable results [5]. 

It is most surprising that the straightforward application of the AIM failed for the 
anharmonic oscillator studied by Barakat [5] since it had been found earlier that the 
approach should be accurate in such cases [2]. 

The main purpose of this paper is to verify whether the AIM gives accurate 
eigenvalues of the perturbed Coulomb model or if its sequences are oscillatory divergent 
as mentioned above. We also discuss the application of perturbation theory to that 
model. 

In Sec. 121 we present the model and discuss useful scaling relations for the potential 
parameters. In Sec. El we apply the AIM to the perturbed Coulomb model directly; 
that is to say we do not convert it into an anharmonic oscillator. In Sec. IHwe outline 
alternative perturbation approaches, and in Sec. El we interpret our results and draw 
conclusions. 
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2. The model 

The problem studied by Barakat [5] is given by the following radial Schrodinger equation 

J2 



H = -7;—, + ^ir:^ + V{r) 



l£_ /(/ + !) 

V{r) = ^gr + Xr'^, (1) 

where Z = 0, 1, 2, . . . is the angular-momentum quantum number, and the boundary 
conditions are ^E'(O) = ^(oo) = 0. We restrict to the case A > in order to have only 
bound states; on the other hand, Z and g can take any finite real value. 
It is most useful to take into account the scaling relations 

E{Z, g, A) = Z'E{1, gZ-\ AZ"^) = \g?/'E{Z\g\-''\ g\g\-\ \\g\-'/^) 

= \'I^E{Z\-'l\g\-^/\l). (2) 

Notice that we can set either Z or A equal to unity without loss of generality, and that, 
for example, E{l,—g,\) = E{—l,g,X). Following Barakat [5] we choose n = 0, 1, . . . 
to be the radial quantum number, and we may define a "principal " quantum number 

u = n + l + l = 1,2,.... 

3. Direct application of the AIM 

Barakat mentions that straightforward application of the AIM does not give reasonable 
results because the sequences oscillate when the number of iteration is greater than 
30 approximately [5]. This conclusion is surprising because it has been shown that 
the AIM yields accurate results for anharmonic oscillators [2], and Barakat converted 
the perturbed Coulomb model into one of them [5] . In this section we apply the AIM 
directly to the original radial Schrodinger equation (jT)). 

By means of the transformation iplr) = (f){r)y{r) we convert the perturbed Coulomb 
model ((1} into a Sturm-Liouville equation for y{r): 

y"ir) = Q{r)y'{r) + R{r)y{r) 
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26' (r) 



Q{r) 



r 



^W = {2[nr)-i^]-^}, (3) 

where 0(r) is arbitrary. It seems reasonable to choose 

0(r) = r'+^e-^^-""' (4) 

that resembles the asymptotic behaviour of the eigenfunction for a harmonic oscillator 
when /3 = or for a Coulomb interaction when a = 0. It leads to 

Q{r) = 4ar - ^^llil + 2f3 

2\ 2 , ,o . m , 2/5(/ + 1) - 2Z 



R{r) = {2X-4a')r' + {2g - Aa(3)r + 

+ 2a{2l + 3) - 2E - /?^ (5) 

We can set the values of the two free parameters a and (3 to obtain the greatest rate of 
convergence of the AIM sequences. From now on we call asymptotic values of a and (3 
to such values of those parameters that remove the terms of R{r) that dominate at large 



r; that is to say: [3 = g/{2a) and a = y A/2. Since the asymptotic values of the free 
parameters do not necessarily lead to the greatest convergence rate [2], in what follows 
we will also look for optimal values of a. 

The Sturm-Liouville equation (jHl) with the functions Q{r) and R{r) ^ is suitable 
for the application of the AIM. We do not show the AIM equations here because 
they have been developed and discussed elsewhere [1,2]. Since the AIM quantization 
condition depends not only on the energy but also on the variable r for non-exactly 
solvable problems, we have to choose a convenient value for the latter [1,2]. Later on we 
will discuss the effect of the value of r on the convergence of the method; for the time 
being we follow Barakat [5] and select the positive root of 0'(r) = 0: 

^° = Va (') 

For concreteness we restrict to Z = A = 1 and n = / = 0, and select 

g = —2,-1,1,2 from Barakat's paper [5]. As expected from earlier calculations on 

anharmonic oscillators [2], the rate of convergence of the AIM depends on the value 
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of a. In order to investigate this point we choose g = —2 because it is the most 
difficult of all the cases considered here. More precisely, we focus on the behaviour 
of the logarithmic error L^ = logl-E''-^-' — E'^^^'^\, where E^'^^ is the AIM energy at 
iteration N and E^"""^^ = -1.1716735847196510437987056 was obtained by means of the 
rapidly converging Riccati-Pade Method (RPM) [9,10] from sequences of determinants 
of dimension D = 2 through D = 22. 

We first consider the asymptotic value a = l/V^. Fig. [T] shows that Ljy decreases 
rapidly with N when A^ <^ 20 and then more slowly but more smoothly for A^ > 20. 
In the transition region about A^ ^ 20 we observe oscillations that can mislead one into 
believing that the AIM starts to diverge. 

Fig. 121 shows that the behaviour of Ln for a nearly optimal value a = 1/2 is similar 
to the previous case, except that the transition takes place at a larger value of A^ and the 
convergence rate is greater. More precisely, Lj^r decreases rapidly with A^ when A^ <^ 50 
approximately as L/^r ^ 0.22 — 0.064A^ — 0.0029A^^ and more slowly and smoothly for 
A^ > 50 as Ltv ~ —6.5 — 0.068A^. Again, the transition region exhibits oscillations. 

TableC]shows the ground-state energies for g = —2, —1, 1, 2 and the corresponding 
nearly optimal values of a. We estimated those eigenvalues from the sequences of AIM 
roots for A^ = 10 through A^ = 80. Notice that the optimal values of a in Table Q 



depend on g and do not agree with the asymptotic value a = y 1/2. Table [T] also shows 
that the AIM eigenvalues agree with those calculated by means of the RPM [9, 10] from 
sequences of determinants of dimension D = 2 through D = 15. 

The rate of convergence also depends on the chosen value of r. The calculation of 
Ljv as a function of ^ = r/vQ shows that Lj\f{C,) exhibits a minimum at C,n and that C,n 
increases with A^ approximately as C,n = 0.435 + 0.005A^ (for g = —2). However, in 
order to keep the application of the AIM as simple as possible we just choose r = tq for 
all the calculations. 



Asymptotic Iteration Method 6 

4. Alternative perturbation approaches 

Barakat [5] first converted the radial Schrodinger equation (PJ into another one for 
an anharmonic oscillator by means of the standard transformations r = u^ and 
$(n) = u^^^'^'^iu^). Finally, he derived the Sturm-Liouville problem 

f"{u) + 2 (^^ - au') fin) + {eu' - 8gu' + 8Z) f{u) = 0, (7) 

where L = 21 + 1/2 and e = ^E — (2L + 5)a, through factorization of the asymptotic 
behaviour of the solution: 

$(m) = M^+ie-°"'/V(M), « = VSA. (8) 

Notice that present a and Barakat's a are not exactly the same but they have a close 
meaning and are related by (y^'as'ymptotic = otBarakat/'^- Since Barakat's application of the 
AIM to Eq. ((7j) did not appear to converge [5] he opted for a perturbation approach 
that consists of rewriting Eq. ((Tj) as 

f'\u) + 2 1^^^ - au^\ f\u) + [eu" + 7 [Sgu^ + 8Z)] f{u) = (9) 

and expanding the solutions in powers of 7: 

00 00 

/(«)=E/^'^(«)7^ e = Ee^'V (10) 

i=o j=o 

The perturbation parameter 7 is set equal to unity at the end of the calculation. The 
series for the energy exhibits considerable convergence rate and consequently Barakat 
obtained quite accurate results with just two to five perturbation corrections [5]. Barakat 
calculated the coefficient e'-'^-* exactly and all the others approximately [5]. 

The model (0) is suitable for several alternative implementations of perturbation 
theory in which we simply write V{r) = Vo(r) + '~fVi{r) and expand the solutions in 
powers of 7. 

If we choose Vo{r) = —Z/r (when Z > 0) and Vi{r) = gr + Ar^ then we can 
calculate all the perturbation coefficients exactly by means of well known algorithms [8]. 

One easily realizes that the perturbation series can be rearranged as 

00 00 
E = Z''J2J2 di,g'\^Z-^^'+^^\ (11) 

j=0 j=0 
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It is well known that this series is asymptotic divergent for all values of the potential 
parameters. 

The other reasonable perturbation split of the potential energy is Vo(r) = Ar^, 
Vi{r) = —Z/r + gr. In this case we can rearrange the series as 

oo oo 

E = Ai/2 J2 J2 Q,ZVA-('+3^)/^ (12) 

i=0 j=0 

One expects that this series has a finite radius of convergence. This is exactly the series 
obtained by Barakat [5] by means of the AIM and, consequently, it is not surprising 
that he derived accurate results from it. In this case one can obtain exact perturbation 
corrections at least for the first two energy coefficients. 

For simplicity we concentrate on the states with n = 0. The eigenf unctions and 
eigenvalues of order zero are 



'01 



n{0) V 2A 



r(/ + 3/2) 



E^T =Y^(2/ + 3) (13) 

respectively. With the unperturbed eigenf unctions one easily obtains the perturbation 
correction of first order to the energy 

p(i)_ Ji + m i\z{2xy/' 

°' (2A)V4r(/ + 3/2) r(/ + 3/2) ^ ' 

that is the term of the series p2|) with i + j = 1. One can easily carry out the 

same calculation for the states with n > using the appropriate eigenfunctions of 

the harmonic oscillator. 

Equation (fT^ yields all the numerical results for e^i in Tables 1-3 of Barakat's 

paper [5]. In particular, E^i = when g = v2AZ/(/ + 1) as in Table 1 of Barakat's 

paper [5]. This particular relationship between the potential parameters also leads to 

exact solutions of the eigenvalue equation (Q). Some of them are given by 



2' ' / + 1' 



v2 \/2XZ 

ES-' =„(2, + 3)-^i^^^j.« = ^, (15) 



where Ni is a normalization constant. 
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5. Conclusions 

We have shown that the AIM converges for the perturbed Coulomb model if the values 
of the free parameters in the factor function that converts the Schrodinger equation into 
a Sturm-Liouville one are not too far from optimal. It is clear that it is not necessary to 
transform the perturbed Coulomb model into an anharmonic oscillator for a successful 
application of the AIM. Our results do not exhibit the oscillatory divergence reported 
by Barakat [5] even when choosing the asymptotic value of a. 

The perturbation approach proposed by Barakat [5] is equivalent to choosing 
the harmonic oscillator as unperturbed or reference Hamiltonian, and if we apply 
perturbation theory to the original radial Schrodinger equation we easily obtain two 
energy coefficients exactly instead of just only one. It is worth mentioning that 
the coefficients calculated by Barakat [5] are quite accurate and, consequently, the 
resulting series provide a suitable approach for the eigenvalues of the perturbed Coulomb 
potential. This application of the AIM to perturbation theory is certainly much more 
practical than the calculation of exact perturbation corrections proposed earlier [7] that 
can certainly be carried out more efficiently by other approaches [8]. 
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Table 1. Energies for some g values calculated by the AIM with TV = 80 and 
RPM [9, 10] with D = 15 



g a 



AIM 



RPM 



-2 0.5 -1.17167358472 

-1 0.3 -0.226186875190871929 

1 0.3 1.33284549226484083 

2 0.5 2.014906226463 



-1.1716735847197 
-0.2261868751908719 
1.3328454922648408349 
2.0149062264617370560 
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Figure 2. Logarithmic error for the energy for g = — 2 as a function of N for an 
almost optimal value a — 1/2 



